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Neutron spin resonance as a probe of superconducting gap anisotropy in partially 
detwinned electron underdoped NaFeo.gssCoo.oisAs 
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We use inelastic neutron scattering (INS) to study the spin excitations in partially detwinned 
NaFeo.gssCoo.oisAs which has coexisting static antiferromagnetic (AF) order and superconductivity 
(Tc = 15 K, Tn = 30 K). In previous INS work on a twinned sample, spin excitations form a 
dispersive sharp resonance near Eri = 3.25 meV and a broad dispersionless mode at Eri = 6 meV 
at the AF ordering wave vector Qaf — Qi — (I 5 O) hs twinned domain Q 2 = (0?1)- For 
partially detwinned NaFeo.gssCoo.oisAs with the static AF order mostly occurring at Qaf = (1?0)? 
we still find a double resonance at both wave vectors with similar intensity. Since Qi = (1?0) 
characterizes the explicit breaking of the spin rotational symmetry associated with the AF order, 
these results indicate that the double resonance cannot be due to the static and fluctuating AF 
orders, but originate from the superconducting gap anisotropy. 


PACS numbers: 74.25.Ha, 74.70.-b, 78.70.Nx 


I. Introduction 

The neutron spin resonance is a collective magnetic ex¬ 
citation observed by inelastic neutron scattering (INS) at 
the antiferromagnetic (AF) ordering wave vector of un¬ 
conventional superconductors below Tc First dis¬ 

covered in the optimally hole-doped YBa 2 Cu 306 +cc fam¬ 
ily of copper oxide superconductors [1] , the mode was also 
found in iron pnictide superconductors at the AF wave 
vector Qaf = (l^O) reciprocal space [Figs. 1(a) and 
1(b)] [SHIS], and is considered one of the hall marks of 
unconventional superconductivity m- Experimentally, 
the neutron spin resonance appears as an enhancement 
of the magnetic spectral weight at an energy Er in the 
superconducting state at the expense of normal state spin 
excitations for energies below it. For iron pnictide super¬ 
conductors with hole and electron Fermi surfaces near 
the F and M points, respectively [Fig. 1(c)] [181 Ej^ 
the mode is generally believed to arise from sign reversed 
quasiparticle excitations between the hole and electron 
Fermi surfaces and occur at an energy below the sum 
of their superconducting gap energies {Er < Ay + Ag) 

[201 ED. 

If the energy of the resonance is associated with the 
superconducting gap energies at the hole and electron 
Fermi surfaces, it should be sensitive to their anisotropy 
on the respective Fermi surfaces [22l [23]. Indeed, re¬ 
cent INS experiments on the NaFei-a^Co^^As family of 
iron pnictide superconductors reveal the presence of 


a dispersive sharp resonance near Eri = 3.25 meV 
and a broad dispersionless mode at Er 2 = 6 meV at 
Qaf = (l^O) ill electron underdoped superconducting 
NaFeo.gssCoo.oisAs with static AF order (Tg = 15 K and 
T/v = 30 K) [T51IT6]. From the electronic phase diagram 
of NaFci-a^Coa^As determined from specific heat [24] . 
scanning tunneling microscopy m, and nuclear mag¬ 
netic resonance (NMR) [26l [27] experiments, we know 
that NaFeo.gssCoo.oisAs is a bulk superconductor with 
microscopically coexisting static AF ordered and super¬ 
conducting phases. For Co-doping near optimal super¬ 
conductivity around x = 0.0175, NaFci-a^Co^^As be¬ 
comes mesoscopically phase separated with static AF 
ordered and paramagnetic superconducting phases m, 
similar to the co-existing cluster spin glass and super¬ 
conducting phases in optimally electron-doped BaFe 2 As 2 
[28l|29]. Since angle resolved photoemission spectroscopy 
(ARPES) experiments on NaEeo.ggsCog.oisAs found a 
large superconducting gap anisotropy in the electron 
Eermi pockets m, the double resonance may result 
from orbital-selective pairing induced superconducting 
gap anisotropy along the electron Eermi surfaces [3T] . 
Upon increasing electron doping to x = 0.045 to form 
superconducting NaEeo.g 55 Coo.o 45 As (Tg = 20 K), the 
superconducting gap anisotropy disappears [32l [33] and 
INS reveals only a single sharp resonance coupled with 
superconductivity [34] . 

Although the superconducting gap anisotropy pro¬ 
vided a possible interpretation [31] , the double resonance 
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in underdoped NaFeo.ggsCoo.oisAs may also be due to 
the coexisting static AF order with superconductivity 
[35l [36]. Since Qaf = Qi = (I 5 O) characterizes the 
explicit breaking of the spin rotational symmetry in the 
AF ordered state of a completely detwinned sample [Fig. 
1(b) and 1(c)] [37] . one should expect magnetic suscepti¬ 
bility anisotropy at Qaf = Qi = (l^O) and Q 2 = (0,1). 
At the AF ordering wave vector Qaf = Qi = (l^O), 
the resonance appears in the longitudinal susceptibility, 
whereas the transverse component displays a spin-wave 
Goldstone mode. At the other momentum Q 2 = (0,1), 
the resonance has both longitudinal and transverse com¬ 
ponents and is isotropic in space. If the resonance shows 
distinct energy scales at Qi and Q 2 , one would expect 
to find a double resonance in a twinned sample as shown 
in Fig. 1(e) [35]|36]. However, one would then expect a 
single resonance of energy Eri at Qi and that of energy 
Er 2 at Q 2 in a completely detwinned superconducting 
sample with static AF order [Fig. 1(f)]. 

To test if this is indeed the case, we have carried out 
INS experiments on uniaxial strain partially detwinned 
NaFeo.ggsCoo.oisAs to study the neutron spin resonance 
at Qi and Q 2 . Instead of E^i at Qi and Er 2 at Q 2 as 
expected from the theory of coexisting static AF order 
with superconductivity [35l |36], we find that both Eri 
and Er 2 are present at Qi and Q 2 as in the twinned 
case. Therefore, the presence of the double resonance 
is not directly associated with the breaking of the spin 
rotational symmetry in detwinned NaFeg.gssCoo.oisAs. 
Instead, our results are consistent with the notion that 
the splitting of the resonance is due to superconducting 
gap anisotropy in the underdoped NaFeo.ggsCoo.oisAs, 
suggesting weak direct coupling between spin waves and 
superconductivity. These results are also consistent with 
polarized neutron scattering data, where the longitudinal 
spin excitations of E^i reveals a clear order-parameter¬ 
like increase below Tc reminiscent of the resonance, while 
the transverse spin excitations of the E^i from the spin- 
wave Goldstone mode have no anomaly across Te m- 

II. Experimental Results and Theoretical 
Calculations 

We prepared single crystals of NaFeo.gssCoo.oisAs by 
the self-flux method [15] and cut a large crystal into the 
rectangular shape along the [1, 0, 0] and [0,1, 0] directions 
(16.11 X 8.41 X 1.31 mm^, ^0.79 g). From NMR measure¬ 
ments m, we know that the tetragonal-to-orthorhombic 
structural transition happens around Tg 40 K, above 
T/v and Tc. Our neutron scattering experiments were 
carried out on the PUMA and BT-7 thermal triple-axis 
spectrometers at the MLZ, TU Miichen, Germany m, 
and NIST center for neutron research (NGNR), Gaithers¬ 
burg, Maryland [38], respectively. In both cases, we used 
vertically and horizontally focused pyrolytic monochro¬ 
mator and analyzer with a fixed final neutron energy of 
Ef = 14.7 meV. The wave vector Q at {qx,qy,qz) in A ^ 
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FIG. 1: (Color online) (a) The electronic phase diagram 
of NaFei-ccCoccAs, where the arrow indicates the Co-doping 
level of our sample. The inset shows the magnetic structure in 
orthorhombic notation, (b) In a 100% detwinned sample, one 
should observe magnetic order at Qaf = Qi = (IT) but not 
at Q 2 = (0,1). (c) The schematic drawings of the Fermi sur¬ 
faces in a detwinned sample, (d) Its twin domain rotated 90° 
away. The arrows mark Fermi surface nesting wave vectors, 
(e) In a twinned sample, one cannot distinguish Qi = (1,0) 
and Q 2 = ( 0 , 1 ), and therefore there should be two resonances 
at both wave vectors, (f) In the AF order and superconduc¬ 
tivity coexisting theory, Eri should appear at Qi and Er 2 at 
Q 2 in a completely detwinned superconducting sample with 
static AF order at Qi. 


is defined as = (ci^aj 211 ^qya 1211 1211 ) recipro¬ 

cal lattice unit (r.l.u.) using the orthorhombic unit cell 
(a ^ 6 ^ 5.589 A and c = 6.980 A at 3 K). In this nota¬ 
tion, the AF Bragg peaks occur at the (1, 0, L) positions 
with L = 0.5,1.5,'' - and there are no magnetic peaks at 
(0,1,1/) [Figs. 1(a) and 1(b)] [39]. We have used a de- 
twinning device similar to that of the previous INS work 
on BaFe 2 -a:NUAs 2 [40]. The samples are aligned in the 
[1,0,0.5] X [0,1,0.5] scattering plane. In this scattering 
geometry, we can probe the static AF order and spin ex¬ 
citations at both Qi and Q 2 , thus allowing a conclusive 
determination of the detwinning ratio and spin excitation 
anisotropy at these wave vectors. Figure 2(a) and 2(b) 
shows the temperature differences in transverse elastic 
scans along the [77,1 — 77, 0.5] and [77,1 — 77, 0.5] direc¬ 
tions, respectively, for NaFeo.ggsGoo.oisAs between 2 K 
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FIG. 2: (Color online) Temperature differences of the elas¬ 
tic scattering below (2 K) and above (41 K) the Neel tem¬ 
perature {Tn = 33 K) at reciprocal positions (a) (0,1,0.5) 
and (b) (1,0,0.5) in a uniaxial strained NaFeo.gssCoo.oisAs 
single crystal. The uniaxial strain along the b axis is ~ 
10 MPa and data was collected on BT-7. By compar¬ 
ing the area ratio of the two peaks, we estimate that ~ 
58% (/(l,0,0.5)-/(0,l,0.5))/(/(l,0,0.5) + /(0,1,0.5))] of 

the crystal is detwinned. The solid lines are Gaussian fits to 
the data. Similar detwinning ratio is also obtained at PUMA. 
After releasing the uniaxial strain, we find that the sample 
returns to the twinned state, (c) Temperature dependence 
of the magnetic order parameters at Qi and Q 2 in unixially 
strained detwinned NaFeo.gssCoo.oisAs. (d) and (e) Constant 
Q-scans at Qi and Q 2 below and above Tc, respectively, in 
the partially detwinned sample. The peaks around ~9 mV 
in the raw data are mostly temperature independent back¬ 
ground scattering, (f) Temperature differences between 5 K 
and 21 K, revealing the superconductivity-induced resonance 
at Qi and Q 2 . The nearly identical resonances at these two 
vectors suggests that the coupling between spin waves and su¬ 
perconductivity is weak in Co-doped NaFeAs. The solid line 
is a guide to the eye. 


and 41 K. By comparing the scattering intensity at these 
two wave vectors, we estimate that the sample is about 
58% detwinned. Figure 2(c) shows the temperature de¬ 
pendence of the magnetic order parameters. Consistent 
with previous data on a twinned sample m, the uniaxial 
strain used to detwinn the sample does not seem to alter 
Tat ^ 30 K and Tc = 15 K. 

In previous INS work on twinned NaFeo.ggsCoo.oisAs, 
superconductivity induces a dispersive sharp resonance 
near E^i = 3.25 meV and a broad dispersionless mode at 
Er 2 = 6 meV at Qi = (1, 0, 0.5) and Q 2 = (0,1, 0.5) [15]. 
To explore what happens in the uniaxial strain detwinned 
NaFeo.gssCoo.oisAs, we carried out constant-Q scans at 
wave vectors Q 2 [Fig. 2(d)] and Qi [Fig. 2(e)] below and 
above Tc. While it is difficult to see the resonance in the 
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FIG. 3: (Color online) Constant-energy scans at the 
low-energy resonance Tri=3.75 meV for detwinned 
NaFeo.gssCoo.oisAs. (a) and (b) The rocking scans across 
Q 2 = (0,1,0.5) and Qi = (1,0,0.5) positions above (21 
K) and below (5 K) Tc, respectively. The wave vector 
independent background scattering increases slightly on 
warming, (c) and (d) Corresponding difference between the 
two temperatures after the background subtraction. The 
solid lines are fits to Gaussians on flat linear backgrounds set 
to zero. 


raw data, the temperature differences between 5 K and 21 
K plotted in Fig. 2(f) reveal a sharp peak at Eri =3.75 
meV and a broad peak at Er 2 = 6 meV, in addition to 
the negative scattering below 3 meV due to a spin gap. 
These data indicate the presence of a superconductivity- 
induced sharp resonance and a broad resonance above 
a spin gap, similar to the results on twinned samples 
m- Surprising, there are no observable differences for 
the resonance at Qi and Q 2 , suggesting that the double 
resonance is not directly associated with the twinning 
state of the sample. 

To confirm the conclusion of Fig. 2, we carried out 
constant-energy scans near Qi and Q 2 at Eri and Er 2 
above and below Tc. Figure 3(a) and 3(b) shows trans¬ 
verse rocking curve scans through Q 2 = (0,1,0.5) and 
Qi = (1,0, 0.5), respectively, at the sharp resonance en¬ 
ergy Erl = 3.75 meV above and below Tc. While there 
is slightly more magnetic scattering at the AF order¬ 
ing wave vector Qi = (1,0, 0.5) compared with that at 
Q 2 = (0,1,0.5) in both the normal and superconduct¬ 
ing states m, the superconductivity induced intensity 
changes shown in Fig. 3(c) and 3(d), defined as the reso¬ 
nance [ 2 ] , at these two wave vectors are indistinguishable 
within the statistics of our measurement. 
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FIG. 4: (Color online) Constant-energy scans through 
the high-energy resonance F^r2=6.5 meV for detwinned 
NaFco.gssCoo.oisAs. (a) and (b) The rocking scans across 
the Q 2 = (0,1,0.5) and Qi = (1,0,0.5) positions above (21 
K) and below (5 K) Tc, respectively, (c) and (d) The corre¬ 
sponding difference between the two temperatures. The solid 
lines are fits to Gaussians on zero backgrounds. 


Figures 4(a-d) summarize wave vector scans at the 
energy of the broad resonance Er 2 = 6 meV around 
Q 2 = (0,1,0.5) and Qi = (1,0,0.5). Similar to data 
at Erl = 3.75 meV, we find that the superconductivity- 
induced intensity gain of the broad resonance is almost 
indistinguishable at these wave vectors, again confirm¬ 
ing the notion that the resonance is not sensitive to the 
twinning state of the system. 

In electron doped NaFei_a,Coa,As, the dominant or¬ 
bital character of the electron pockets would be d^yj^z 
at (1,0) and dy^jxy at (0,1) in the Brillouin zone [Fig. 
1(c)] [3a|32l|33]. The orbital character of the hole pocket 
is dxzjyz' If fhe superconducting pairing amplitudes are 
highly orbital dependent,z.e., l\xy 7^ ^xzjyz^ fhe super¬ 
conducting gap can be anisotropic along the electron 
pocket and this gap anisotropy gives rise to a splitting of 
the resonance peak m- Such an orbital-selective pair¬ 
ing scenario is consistent with both ARPES measure¬ 
ments m and INS results in twinned samples US]. 

In the uniaxial strain detwinned sample, the degen¬ 
eracy of the dxz and dyz orbitals is lifted and, corre¬ 
spondingly, the Fermi surface is distorted [30]. To in¬ 
vestigate whether the double resonances in the orbital- 
selective pairing scenario still exist in the presence of a 
splitting between the dxz and dyz orbitals, we calculated 



FIG. 5: Imaginary part of the dynamic spin susceptibility x" 
in the superconducting state of a multiorbital t —Ji — J 2 model 
with a nonzero splitting between the dxz and dyz orbitals, 
e = 0.02 eV. Two resonance peaks are present in each of the 
Qi and Q 2 wave vectors. In the calculation, J 1 /J 2 = 0.1 is 
taken such that the pairing amplitudes show strong orbital 
selectivity. See Ref. m for details of the model and the 
method. 

the imaginary part of the spin susceptibility x" in the 
superconducting state from a multiorbital t — Ji — J 2 
model with an orbital splitting term e = Uxz — "^yz ED- 
Our result for a strong orbital selectivity is presented 
in Fig. 5. We find two resonance peaks at each of the 
wave vectors Qi and Q 2 for a nonzero splitting e. The 
intensities of the counterpart peaks at Qi and Q 2 are 
comparable. At each resonance peak, there is a relative 
shift of the resonance energy between the Qi and Q 2 
resonances. This shift is proportional to the splitting e. 
The calculated double-resonances feature at both Qi and 
Q 2 is qualitatively consistent with the experimental ob¬ 
servation in the detwinned sample. The experiment can 
not resolve a relative shift of the resonance energy. This 
could be because either the splitting e is small in the de¬ 
twinned underdoped compound, or the coupling between 
the superconductivity and the splitting e is rather weak. 
Further comparison between theory and experiments is 
needed to fully settle the issue. 

III. Conclusion 

In conclusion, our INS experiments on partially de¬ 
twinned NaFeo.QssCoo.oisAs reveal the presence of two 
resonances at each of the wave vectors Qaf = Qi = (I 5 0) 
and Q 2 = (0,1). This is different from the scenario where 
the two resonances are due to the coexisting AF order 
with superconductivity [35l |36]. Instead, the data are 
qualitatively consistent with the proposal that the dou¬ 
ble resonances originate from an orbital dependence of 
the superconducting pairing. Our results provide further 
evidence that orbital selectivity plays an important role 
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in understanding not only the normal state but also the 
superconducting pairing of the multiorbital electrons in 
the iron pnictides. 
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